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ABSTRACT 


The classic single echelon inventory model is restruc- 
tured as a two echelon problem in which demand and resupply 
are deterministic. Using cost minimization as the objective, 
three models are developed which address the problems of 
(1) no stockouts allowed (EOQ), (2) backorders allowed, and 
(3) finite production with no stockouts allowed. General 
solutions for the optimal policy are obtained in the EOQ and 
finite production models. In the backorder model, the ana- 
lytical argument is limited to the case in which only time 
dependent backorder costs occur. Algorithms are developed 
for solving problems for all three models, and selected para- 


meter values are used to test the behavior of the models. 
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IL... ШЫТКОБ ШАҚЫ 


Traditional inventory analysis addresses the single 
echelon inventory problem. In this problem, customer demands 
are received and issues are made from a single outlet. The 
outlet in turn is replenished from a single source. 

Unfortunately, most real world inventory systems are not 
this simple. A major manufacturer, for example, usually has 
a production and distribution system that includes (1) fac- 
tories, (2) factory warehouses, (3) regional warehouses, 

(4) local warehouses, and (5) retail outlets [Hadley and 
Whitin]. Such a system is an example of a multiechelon inven- 
tory system. Each level of the system is a separate, distinct 
echelon. Figure 1 depicts graphically this organization. The 
United States Navy with its organization of inventory control 
points, supply centers, supply depots, shop stores, tenders and 
underway replenishment ships provides an example of a larce 
multiechelon inventory system. 

In defending their single echelon models Hadley and Whitin 
argue that, even though most real world systems are multi- 
echelon, it is often true that the system need not be treated 
as multiechelon. They contend that a different organization 
frequently operates each level (echelon) of the system. In 
the example of Figure 1, a factory and its associated ware- 
houses might be operated by an equipment manufacturer whereas 


the regional warehouses might be separate privately owned 
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Figure l. Typical Multiechelon Inventory System 


distributorships not under the administrative control of the 
manufacturer. In like manner, the local outlets might be 
individually owned retail stores. Thus we would have a sys- 
tem comprised of several single echelons which are linked by 
a physical dependence but separated by administrative 
controls. 

While it is clear that there are many examples of multi- 
echelon systems that can be treated as a series of independent 
single echelon problems, it should be equally clear that there 
is a large class of systems which must be treated as strongly 
dependent multiechelon. 

It is the objective of this thesis to develop and investi- 
gate the behavior of multiechelon models constructed from the 
classic single echelon inventory models. Specifically, the 
following deterministic single-echelon models of Hadley and 


Whitin will be restructured as two-echelon models: 


> 





a. No stockouts allowed (Economic Order Quantity), 

b. Backorders allowed, 

c. Finite production with no stockouts allowed. 
Like their single echelon counterparts, these two-echelon 


models will be single item, single source models. 





ТТ. -THE TWO ECHELON MODEL WITHANO STOCFOULS 


In any study of inventory theory, one of the first models 
investigated is the deterministic demand model with no stock- 
outs allowed; often referred to as the Economic Order Quantity 
(EOQ) model. While it is true that complete deterministic 
demand is almost never known, it is felt that the mathematics 
of this model provide a good starting point for any inventory 
analysis. Further, it is felt that the deterministic model 
will provide an insight into the operation of a stochastic 
demand model. It is considered appropriate, therefore, to 
begin this multiechelon analysis by considering a two echelon 


extension of the classic EOQ model. 


A. MODEL FORMULATION 

The two echelon model is based on the following 
assumptions: 

1. The upper echelon always replenishes its supply from 
the same outside scurce. The lower echelon always replenishes 
its supply from the upper echelon. 

2. The upper echelon is always able to meet the demand 
of the lower echelon. 

3. External customer demand always occurs at the lower 
echelon. 

4. The external demand rate is deterministic, continuous, 


and constant with a value of A units per year. 





5. Neither backorders nor lost sales are allowed. 

6. Procurement lead time (PLT) is negligible. 

7. The order quantity at the upper echelon, Q2, is an 
integer multiple of the lower echelon order quantity, Ql. 

Let n = Q5/04 ; then this assumption requires n > 1 and 
integer. 

The optimal policy throughout this thesis will be that 
which minimizes average annual variable system cost, subject 
to the constraint(s) of the model. The form of the average 
annual variable system cost will be developed by first deter- 


mining the system cost per cycle, and then dividing this cost 


by the cycle length. The total system variable cost per cycle 


is the sum of the individual echelon variable cycle costs, 

K} (c) and Ко (с). The cycle length is defined as the time be- 
tween receipt of two successive orders at the upper echelon 
(echelon two). 

The cost per cycle at echelon two is comprised of an 
етапа cost, an inventory holding cost (IHC) and a purchase 
cost. The ordering cost is assumed to be a constant cost per 
Exder which includes the administrative costs associated with 
inventory review and order (contract) preparation. This cost 
is independent of the quantity on hand or on order. The IHC 
includes a warehousing cost, an obsolescence cost and a fore- 
gone opportunity cost. It is assumed to be a €unction of the 


inventory on hand and to be expressable by 


T 
IHC = 2с а (ван, 


TO 
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Figure 2.  Two-Echelon EOQ Model. 


where І, is the average cost per dollar invested in inventory 
per unit time, C is the item unit cost, T is the cycle length, 
and ef, (t) is the on hand inventory at echelon two at time t. 

The purchase cost is assumed to be independent of the quantity 
produced, and can be expressed as QC where "p is the quantity 


procured. 


The variable cost per cycle at echelon two can now be 


expressed as 


HE 





К. (Е) ues НА гас Ја 008 (2-1) 
The height of the on-hand inventory at echelon two is 
Q, =- Q) = мо) 
at the beginning of the cycle since 0, is ordered, arrives 
immediately, and is partially used to meet an immediate demand 
of Q1 units. The next demand at echelon two is for an amount 
о, which occurs Q4/A time units later. Similarly, another 
demand of Q1 Occurs after Q/A additional time units have 
passed. The on-hand level at echelon two continues to de- 
crease in steps of Qi until it reaches zero. It remains at 
zero for the last Q4/A time units of T. The area under the 
on-hand inventory curve for one cycle of echelon two is 
therefore 
nl 


E Q 
= 1 „enin=1]..22 
St. et E NN EE 


m=1 


From (2-1), the total variable cost per cycle ас echelon 


two is 
2 
І.С n[n-1] Q1 | 
К. (с) = A, + CQ, + 55 (2-2) 
From the assumption that 0, = nO, , where n is integer, 


it follows that the cycle length at echelon two is equal ton 
reorders. at echelon one. Therefore, from Figure 2, 


Q nQ 
те са. (Bas) 


The average annual cost at echelon two is obtained by 


dividing (2-2) by (2-3) which yields 





AJA 1,C[n-1105 
K = —— + Cl + —— , 
2 Qj 20, 


ШЕ, as a’ funetion of 01 апа п, 


АА ICQ, [n-1]}] 
_ 2 2 = 
Ko == nO, T CA Е сео еи • (2-4) 


Development of the cost per cycle at echelon one follows 
an argument that is analogous to the echelon two development, 
except that there is no purchase cost incurred at echelon one. 
Also, from the assumption that о. is an integer multiple of 
О, , it follows that in each cycle there are n reorders at 
echelon one, i.e., к. (с) is linear in n. Thus the cost per 
cycle can be written as 

пт. со: 
Hc ы s (25 
When (2-5) is divided by (2-3) the average annual variable 


most 15 


А.А I CO 
Е 
" = 5 . (2-6) 





The total average annual cost of the system, K, is the sum 
of (2-4) and (2-6). 


A4À (O ALA I;C[n-1]Q 
4 el + Een CA + = I 
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ЈЕЛ 


The model has assumed that the unit price is independent 
of the quantity ordered. Therefore the C term in (2-7) is a 
constant, and can be dropped from (2-7) resulting in (2-8) 


being the average annual variable cost. 


ЕЗ 
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i 
ЕН OPTIMAL POLICY 
Determining the optimal inventory policy E finding 
the values of Q1 and n (call them QF and n*) which minimize 
(2-8) for a given set of model parameters. To find this 
policy, rewrite (2-8) as 


A 


[I, + I.,[{n-1]]co 
K (Q4 ,n) EP ГАЈ + =) i i. 


À 

=— + ---- (2-9) 
Q1 2 

Then since K (Q4 ,n) is convex in 04 for fixed n, QF (n) must 


satisfy the equation 


JK 1. за. А [Ii + 1, [1-1]]С " 
90, 1 n m 2 
Or 
А, 1/2 
+ 2A [А] + a 
Q1 (n) = ТЕ + ЕВЕ ИЕ ‘ (2-10) 


Ester substituting (2-10) for Q4 in, (2-9) 2ondzeolleeting 


terms, (2-9) takes the following form: 


A, Ша 
K (Q* (n) ,n) = {2A [A, а то (е e ; 
Or 
2ХА Т.С 
K (03 (n) ,n) = БЕС Foe 
ҖЫ 1/2 
+ 21 A I,C[n-1] + 1) . (2-11) 


Note that the first term in (2-11) is equivalent to ке 


ИВ. (2-11) [Hadley and Whitin]). Іп fact, it is readily scen 


that each term in (2-11) contains a form of this classic 


кола bet 


т 
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Then (2-11) can be rewritten 


к2. 2 
K(n) = о, ‚п) = ша + таа + K],[n-1] + K [—— 


"EN ) 


It is immediately seen that (2-12) is positive for all 


feasible values of n (n > 1). Further, 
m Kn) = [к^ + 3122 , 
nob С 


and for very large values of n, 


Ва 





1% K(n) > Ky > n > 0 
K 
oK (n) 12 
on 24172 
2 -K 
3. д En E 11%; «meo 
n 4 „3/2 


If the assumption that n is an integer is temporarily re- 


laxed the slope of K(n) can be readily investigated by 








considering 
2 2 
mor. 52 2 
LK 
2 2 12 
К n n 
n ^ "c к. 
2K (n) 
If (2-13) is positive when n = 1l, then it follows that it 
is positive over all n, and the optimal value of n is n* - 1. 
If (2-13) is negative when n = 1, then there exist values of 
n such that K(n) is less than K(l). However, since it has 


been shown that (2-13) is positive for very large n, it follows 


that n* is finite, and satisfies 


1 





on sr T 
Or 
2. „2 1022 
ne (421. 22 > 22 | (2127) 
nm. 
A further consideration of the case where (2-13) is nega- 
tive at n = l reveals that 
2 
8.5 > 0 
on 


over the range of n values from 1 to beyond n*. Let the value 


of n where 


be denoted by n. Then in the range l <n < n the fünctlon 
K(n) is convex. For values of n > n, (2-12) is not convex. 
However, (2-13) will alvays be positive. Therefore n*, given 
by (2-14), is the optimal non-integer value of n. The reader 
will note that the above arguments indicate that (2-12) is 
pseudo-convex. 

Since K(n) is pseudo-convex over all n, it follows that 
K(n) is also pseudo-convex on the integer values of n. If n* 
is not an integer, then the optimal integer value of n will be 
either n (the smallest value greater than n*) or n (the 
largest value less than n*). The integer optimal will be the 
value of n corresponding to К* = min [K(n) , K(n)] . 
In summary, the optimal inventory policy ts found as 


follows: 


16 





2 25 2 
T3 Compute Кој! Koos Ку, 
2 Test КЕ = к“ > На If this inequality. holds then 
Í 21 22.750532 


n optimal = 1. Go to step 8. 


J: Та кең = Kee < KT, then “compute n* ивїпс (2-14), 

4. If n* is an integer then n optimal - n*. Go to 
step 8. 

5. If n* is not an integer, compute n and n. 


6. Compute K(n) and K(n) from (2-12). 


7. K* = min (K(n,K(n)) and n optimal is the value of n 


corresponding to K*. 
8. Compute 03 from (2-10). 
C. NUMERICAL EXAMPLE 


AS an illustration of the algorithm, consider the follow- 


ing problem where 


I, = 0.75 A, = $200.00 
I, = 0.50 А = 100 units/year 
А, 925 100 С = $100/unit 
Begin by computing кї, =: 55 

ШІ 
в ос 
K 

Next compare A “o кан eode x 109 > ке, - 2099 107, 


ЦЕ] 





so compute n* using (2-14). The value is 2.0 whieh is 
integer. Therefore n optimal = 2. Using (2-10) the optimal 


order quantities are 


Q1 = 166.67 units 
and 
05 = 166.67 x 2 = 333.34 units. 


D. PARAMETER ANALYSIS 

The primary objective of this analysis was to observe the 
effect of parameter variations on optimal n. The values of 
four of the five parameters were fixed and the value of the 
remaining parameter was allowed to vary. The process was then 
repeated for each of the other parameters so that the effect 
of each could be observed. 

Representative graphs of optimal n versus a given para- 
meter are illustrated in Figures 7 through 10. In Figure 7 
it is seen that n optimal is inversely related to А}. How- 
ever, as noted in Figure 8, there is a direct relationship 
between n optimal and A,. Further, a comparison of Figures 
7 and 8 reveals that, for the parameters selected, n optimal 
is much more sensitive to A, than toO Ај. 

Figures 9 and 10 indicate the relationship of n optimal 
ЕО Ii and Io, respectively. The results indicate that n 
optimal is directly related to Ii, inversely related to 1, 
and that n optimal is more sensitive to a change in holding 


cost at echelon two than at echelon one. 


18 





An analysis of the effect of demand rate on n optimal 
revealed that n optimal was completely insensitive to A over 


the range 5 € À < 10,000. 


13 





III. -FINITE РБКОРИСТТОММАТІЕНЕРЕБЕФТИОИ РБИЕ БЕ 


An additional two echelon variation to the deterministic 
models treated [Hadley and Whitin] addresses the problem of 
finite production. This model would have application in any 
situation where the upper echelon manufactures as well as 


warehouses a given product. 


A. MODEL FORMULATION 

The model, shown graphically in Figure 3, assumes that 
echelon two produces the needed material at a constant con- 
tinuous rate of Y units per unit time (W >A). The cost A, 15 
now considered to be a set up cost which is incurred each time 
a new production eat: initiated. This set up cost is inde- 
pendent of the quantity produced. The rest of the assumptions 
are identical to those of the EOQ model. 

The model will, as in Section II, seek optimal values of 
Q4 and O, which minimize the average annual system cost. This 


Cost is the sum 


K = K, +k. , 


where к. is the average annual cost at the ith echelon. Іп 


Section II it was shown that 


A 
Lp 
Ki = O T 51 СО 5 


To develop к. it is necessary to define the cycle length 


as 


E 
|| 


ty 4 Ще + Ta ; | | 2 (3-1) 
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where ty is the time from the beginning of the cycle until the 
start of production; Тр is the production length, and Та 15 
the time required to deplete the inventory at echelon two 


after production stops (see Figure 3). Note that 


Ex EN 
Да 27 Ды 
ША 2 22 
р А О ' 


апа 


A 





Define i as the number of times that the quantity Q4 is 


demanded at echelon two during the production period, Tp: TE 


can be seen that 


i = (Y) (3-2) 


{ Qi 
where t = Қаны 


\ 


As in Section II, the variable cost per cycle at echelon 


two, K,(C), will have the form 


t. +T FT 
l p d 
КЕ о PAS dt (33) 
1 


From Figure 4 it follows that 


и 
р а | 
с (t) at = Е(1) + g(i + 1) + h(i + 2), 


1 


ша 





where f(i) is the area under the first i saw teeth, g(i + 1) 
is the shaded area shown in Figure 4 and h(i + 2) is the area 


remaining under the curve. 
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Figure 4. Finite Production Model, Area Under the Inventory 
Curve at Echelon Two. 


From Figure 4 it can be shown that 


02 11-11 07) . от [4-Х 1 [1-1] [1-2] 


fi) = == + (3-4) 
ev 2 222 
апа 
2 | 
| Q1 [n-i] [n-i-1]} 
h(i) = Se (3-5) 
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Figure 5. Finite Production Model. Area Under 
1 + ыг Saw Tooth. 
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Substituting (3-4), (3-5) and (3-6) into (3-3) yields 





Tacos .с01913-11 
КІСІ ы Ф221 


2 2) 2) 2 


1,C0* [n-i] [n-i-1] | 1,C0* [n-1]* 


24 24 


Deu DER 
_ 1С01У11-1) 4 T,cQ, liln-i]] 


2) 2 A 


IjC[i-1] [j3A] 
—€———— , (3-7) 
À 


After dividing (3-7) by the total cycle length and collecting 
terms K, can be written 
А.А 1,00, [1-1] I4CQ4A [n-2] 


Ko = no, F жеу = I ен ° (3-8) 


The system average annual variable cost is found by summing 


K, and K 


1 2. 


А 


А [1,+1, [n-1] - 1,7 


2 
2 


[n-2]}co, 
K (Q4 ,n) = [A +] + - 


(3229 
ES OPTIMAL POLICY 

The development of an algorithm for determining the optimal 
inventory policy in the finite production model follows an 
argument that is identical to the one discussed in Section II. 
Since (3-9) is convex in Qi for fixed n, ОТ (а) must satisfy 


OK (Qn) 


30, S 
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Or 
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2 [A 4-5] 1/2 
Оя (п). = онаст к | (3-10) 
C [Ii*I5[n-1] = I,gin-21 
when (3-10) is substituted for 0, in (3-9) then 
42 
HN Е + = + к“, e 18-2] 
K2,l[n-1] - Ма-211|5/2 
+ Se 


where к^. о. Со 
2 18 2) 


It is apparent that (3-11) is positive for all values of n. 
Further, 


Е 2 D DER 
lim K(n) е + K + KT — 
n>] TI | 247 120 22106. 


and for very large values of n, 
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If the assumption that n is an integer is temporarily 


relaxed the slope of K(n) can be investigated using 
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ТЕ (3-12) is positive when n = 1, it is also positive over 
all n and it follows that the optimal value of n is n* - 1. 
If (3-12) is negative when n = 1, then there exist values ene 
n such that K(n) is less than K(l). However, Since it has 
been shown that (3-12) is positive for very large n, n* must 


be finite and satisfy 


ок 








п 
Or 
n* = Р (359) 
where 
^2 М 2 -2Х 
Еа ЛЕН! 
and 
202 ЕЕ 2 ¿YA 
meine: 
A further consideration of the case where (3-12) is nega- 
tive for n = 1 will show that 
2 
9 a > 
on 


over the range of n values from l to beyond n*. Let the value 


of n where 


be denoted by n, and it follows that in the range 1 <n < n 


Kin) is convex. For values of n ^n K(n) is not cener TERE 


27 





(3-12) is positive. Therefore n* is the optimal non integer 
value of n and (3-11) is pseudo-convex. 

Since (3-11) is pseudo-convex over all n, it follows that 
it is also pseudo-convex on the integer values of n. So if 
n* is not an integer, then the optimal value of n will cor- 


respond to 
K* = min[K(n), K(n)]. 


The optimal policy algorithm can now be written. 


and Ke 


2 ^2 
l. Compute Кој! Koos 12 


2 22 ^2 не. 
3 Tf Ke > 52 4- R2 then compute n% usin (3-13) 
: 21 22 12 р J ; 


4. If n* is an integer then n optimal = n*. Go to step 


5. If n* is not an integer compute n and n. 
6. Compute K(n) and K(n) from (3-11). 


7. K* » min[K(n), K(n)] and n optimal is the value of n 
eor responding to K*. 


8. Compute QF Dom ODE 


6. SPECIAL CASES OF y 





When Y =A 
__АјА_ 1,00) А 1,CQy 
о E 
ji = 
and 
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A [I.4I.]C0 
jid C = + ж м (3-14) 
: | 


Dico 


Thus it is seen that if the production rate is equal to the 
demand rate, the optimal policy is to start the production 
line and never let it stop. 

As one would expect, when у > o 


Pun 1,00. АДА 





1 2 
Tim Kn) =- =} + — 
y> Q1 2 пој 

I CQ} [n-1] 
ШЕ" 


and it is immediately seen that this model degenerates to 


the EOQ model of Section II. 


D. NUMERICAL EXAMPLE 


Consider the problem where 


т = 0350 er = SNO dC 
I, = 0.20 А, = ТЫРЫП С 
А = 25 units/year А, = 2200/83 t 
ЕЕ ЕЛЕ 
2 2 да 
Hirst Compute Rig! Ко? апа Ко“ 
2 E 5 2 T 5 
Ko, = 30 x ПО Ка = 1252210 
2 _ 5 22 E 4 
Коу = 2.0 x 10 K55 = 2227 52 110 
2 а 4 2 кй. 4 
КІ: = 5.0 х 10 Ко = JOE 
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Next compare Ka with [os ү Ка, К = 5.0 x 102 ? Kas + = 1.0х io 
so compute n* using (3-13). The value is 3.64 which is not an 
integer. Therefore n = 4 and n = 3. 
Prom Pe 
К(п) = $685.00, 
К(п) = $694.00, 


and n optimal is 4. Using (3-10) 


Qi = 7.2 Unies 
апа 
Q5 = 28.8 units 


E. PARAMETER ANALYSIS 

Figures 1l and 12 illustrate the effect of production rate 
on optimal n. If the holding cost at echelon two is small 
then optimal n was relatively insensitive to a change in 
production rate (Figure 11). However, the results of Figure 
12 show that this sensitivity increases as I, increases. 

The responses of this model to parameter variations were 


identical to the responses observed in Section II. 
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IV. THE TWO ECHELON MODEL WITH BACKORDERS PERMETTED 


A. MODEL FORMULATION 

The EOQ and finite production models were based on the 
assumption that all demands would be immediately satisfied. 
This assumption will now be relaxed in order to investigate 
a System in which all demands are ultimately satisfied, but 
where it is permissible to accumulate backorders at the lower 
echelon. No backorders are allowed in echelon two. 

In this model demands which occur when the lower echelon 
is out of stock are backordered against future procurement. 
When the procurement arrives these backordered demands are 
met, and any excess quantity is placed in stock. Each back- 
order at echelon one results ina cost of the form m + п, 
where m is a fixed charge and " is a time dependent cost. 

As noted, if there were no costs associated with incurring 
backorders [Hadley and Whitin], then it would be optimal to 
never have any inventory on hand. Conversely, if the back- 
orders are sufficiently expensive, then the optimal policy 
would be to never incur any. However, for an intermediate 
range of backorder costs, it will usually be optimal to incur 
some backorders toward the end of the cycle. 

With the exceptions noted in the preceding paragraphs, the 
model is predicated on the same assumptions used in Section II. 

Define s as the number of backorders at echelon one when 


nnoorder ie received... Then from Figure 6 it ftollewe trhat (he 


oy 





on hand inventory at echelon one immediately after receipt of 
an order is (Qi - s) units. It is assumed that Q4 > в. Then 


the echelon one on hand inventory varies from (Qi -'5) to 


zero. 
Echelon II 
До ч) 
е; 
Y 
fi 
4) 
> 
я 
A 
го 
с 
© 
“С 
с | 
О Time 
M s LM 
À 
> 
$4 
О 
43 
E 
0) 
3 Echelon I 
B 
rj 01-5) 
с 
© 
с 
С * 
с Time 


о 
A 


Figure 6. Two-Echelon Backorders Allowed Model. 


The optimal policy continues to be that which minimizes 
average annual variable system cost subject to the model 
constraints. Since the model does not permit backorders at 
the second echelon, the average annual variable cost at 
echelon two will be identical to equation (2-4), which for 
convenience is restated here. 
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ДЕ) Т „СО. [n-1] 
E 2 2” _ 


Development of the echelon one variable costs follow 


identically the single echelon backorders permitted model of 


Exdley and Whitin. In the interest of brevity, this develop- 


ment will be omitted and the results (equation 2-17) [Hadley 


and Whitin] merely stated. 


is 


2 
А» Т.С[О--5] ^ 2 
mE D е 4: TS к 
Ki = 23 + по + Qr (ТА 5 Tem ] (4-2) 


Summing (4-1) and (4-2), the average annual variable cost 


2 
АД АЛ I;,CÍQ,-s] 032 
ill 1 1 1 TS 
К = = + — + - + — [ns + —— (4-3) 
Q1 no, 20, Q1 2 
1200. [1-]] 
TE oe i 


OPTIMAL POLICY 


The optimal solution seeks the values of the decision 


variables QT, S*, and n* which minimize (4-3) subject to the 


constraint that all demands are ultimately satisfied. To do 


this, note that (4-3) can be rewritten as 





LEO 2 
AA 1 TÀS Х 5 
К Ap И ер : (4-4) 
Q1 2 1 01 201 
where 
X - ICs 7, 
A 
> 2 
А = Ау +=, 
апа 
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[nz її] 


it follows that 


S = 
or 

om = 
Next, 

ASAS 

904 


BESHults in 


Qi = 





PO TA 
ТА + XS 
Е (4-6) 
ТС 
-2M TE 2 TAS XS _ 
DEM EE 2e вот РАО 
204 204 201 
5 172 
РЕ Ре + XS (4-7) 
IC 


To solve for 07 and s* in terms of n and the system parameters, 


substitute (4- 


7) into (4-5) and (4-6). After collecting terms, 


the resulting equations are 


QT (n) 
and 


веки) 


1/2 
Е 2A X - Tr 
MES A (4-8) 
NS IMS 
1 
1/2 
ek tone = E тА T 
В a ml c Ја 
то пс 


The expressions (4-8) and (4-9) can now be back substituted 


Tor Q4 and s in (4-4) to yield an expression for the average 


* 
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annual variable cost, K, as a function of n and the system 
parameters. After regrouping the terms this expression can 


Be written 


гр 12 а ос 
[IC[1,C+]-1,C ] [2A A[I, Cr] 77 А“) I Cm 
Ee DE — ——— — tT + 
п.с + T] JT 


(4-10) 

Equations (4-8), (4-9), and (4-10) are valid only if 
2AMIC - ns 2 0. Otherwise an analysis similar to that of 
Hadley and Whitin is required for the case when m and T are 
positive. Such an analysis is beyond the scope of this thesis. 

Because of the complex form of K(n) given by (4-10) as a 
function of n no attempt was made to evaluate whether or not 
K(n) was convex or pseudo-convex in n when т апа T are both 
positive. However, an example is presented later in which an 


optimal value of n > 1 is obtained for T and Т positive. 


If пл = 0 (4-10) reduces to 








1/2 
IC 
К(п) = [2AAIC - 2MI,C =! , (4-11) 
To expand (4-11), let 
ке, = 2А.Т.С (6) 
1j E 
Then, 
k X 
=. oe 2 2 _ 21 
К(п) = -172 Ky 1x + Кох [п 1] + = 
KS 4X [n-1] 2 Ko. poe a 
y E АКО. Се е 
n Jun EN E n í 
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and since x = І.С + Т, K(n) reduces to 





1 
2 ~ 2 ~ 
K ТІ K T 
ка) = |== + SL + Kf [n1] 
1,C+T n [1,C+T1] 
к 1/2 
+ = z (4-13) 


n 


Following the analysis of Section II, it is immediately 
seen that (4-13) is positive for all feasible values of n, and 
that 


lim K(n) = 
nel I.C + 


Also, as n gets very large 


1/2 


1 K(n)>»K,, n О 
К 
ок LZ 
2. 7 > 0, 
on 2112 
Кол р. 
а 2n 7? 


Assume, for the sake of argument, that n is continuous 


Or all n > l. Then 


2 2 
en 2, O 
—————— € Ki. + — 
2.2 12 2 
2 [I.C-tT]n n 
ð K E 1. 
—— Se cal) (uS (4-14) 
дп K (n) 
and it is obvious that for n = 1 (4-14) is non negative if and 
omy if 
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| TC 
2 2 2 1 


If (4-15) holds, them (4-14) з5 positive for alin. tars 
implies that over all n, the slope of K(n) is positive, and 
m* = 4. 

IF (4-14) is negative at n = 1, then there exist values 
of n such that K(n) « K(1). However, it has been shovn that 
for very large n, (4-14) is positive, and it follows that n* 


15 finite and satisfies 


ТТ 
Or 
N 272 
каут - КА (11С + T] 
n* = ee — (4-16) 
К.С + m] 
To complete the argument, note that when (4-14) is 
negative at n = 1, 
2 
25>0 
on 


over the range of n values from 1 to beyond n*. Let the value 


of n where 


be noted by n, and it follows that K(n) is convex over the 
range 1 сп « n. For values of n » n, (4-13) is not convex, 
but in this range (4-14) is always positive. Therefore n* 


given by (4-16) is the optimal non integer value of n. The 
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reader should again note that the above arguments indicate 
that (4-13) is pseudo-convex. 

Since K(n) is pseudo-convex over all n, it is also pseudo- 
convex on the integer values of n. Therefore, if (4-16) does 
not yield an integer value then the optimal value of n can be 
found by evaluating K(n) and K(n) as in Section II. 


An algorithm for finding the optimal policy (in the case 





where T = 0) can now be stated. 
2 2 2 
1. Compute Кој! Кор! апа Ку : 
2 2 2 e 
2. Test Кот Са [Ky 5 + K5] [1 + ^ ] . If the inequality 


holds then n optimal is 1. Go to step 8. 


2 2 2 m | 
ce. If Ко > Ку, + Ко.) (1 + кш then compute n* using 


4. If n* is an integer then n optimal = n*. Со to step 8. 

5. If n* is not an integer compute n and n. 

6. Compute K(n) and K(n) from (4-13). 

7. K* = min[K(n) and K(n)] and n optimal is the value of 
n corresponding to K*. 


8. Compute ОТ from 


> 
24 [A4 + =a 1/2 
QT (n) = 


C[1,1,C[n-1] + Т.т + 1,1[n-1]] 


1 


ana s* from 


u * 
IC 01 1/2 


1,C + T 


Оа 
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C. NUMERICAL EXAMPLE 


Consider the problem where 


т = 022745 m - 055000 

I, = 0.25 " - $100.00/year 
A, = $25.00 А = 100 units/year 
A, = $100.00 C = $100/unit 


| | 2 Е 
Begin by computing Ki = e e 


2 к 6 

2 Y 5 
Ко? = 5.0 x 10 

2 _ 5 
Кт, = 725 x TO 


dec 


2 2 2 1 Р 
Next compare K5i < ІК1- Е Koo] [1 + —]. The inequality does 


T 
not hold so compute n* using (4-17). The value is 1.69 which 


is not an integer. From (4-13) 


K (n) $735.00 


and 


K (n) $620.00 


Therefore optimal n 1. The values of От ‚ О®, апа s* 


2 
are now found to be 


0% = 12.9 units 
05 = 12.9 units 
5% = 8,16 units. 
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D. PARAMETER ANALYSIS 

As noted previously, computational difficulties preclude 
an analytical treatment of the general form of the, backorder 
model. However, the response of the model to parameter 
variations can be approximated. This can be accomplished by 
considering separately the two functions (3-1) and (3-2) and 
ignoring, for now, the integer restriction on n. 

Equation (4-1) can be rewritten 


2 А.А А IjC[O, - Q,] 
= ња ЕЕ: 
2 Q, 2 
and since the second term is never negative, it is obvious 
ehnat K, is convex in e . Using the Hessian it can be shown 


that (4-2) is convex in Qi and s if and only if 


2A,A [1,C + 7] > [mé . (4-18) 
Then 
АЈА ТС [© - EE 1 най 
K == ТТ + mE т T О [745 + ——] 
l 1 1 


Ж ТЕ [О = а 
2 2 2 č 1 (4-19) 
Q^ 2 


+ 
is convex if and only if (4-18) holds. 
It can be shown that when (4-19) is not convex the optimal 
inventory policy will be realized when backorders are not 
allowed (s* = 0). Thus if (4-18) does not hold, the problem 


can be solved by setting s* - 0 and using the algorithm 


described in Section II to get Or and n™. 
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For the case where (4-19) is convex, optimal 014, Q^ and 


S can be found from 


= he — = Оха (4-20) 


Equation (4-20) will be satisfied when 


2. [nÀ, * A] [IC € 5) - nmr“ v2 

0% = Se) LL. ee: 5 = , (4-21) 
n[I,c кт + 1112С [n-1] + Тост {п=1]| 

1 СОТ (а) - ТА 
5* (п) ч -----------, (4-22) 
т. С + TI 
and 
Q5 (n) = поў (п) , (4-23) 


and the optimal policy can be found from the following 


algorithm. 
I Set n=- l- 


2. Compute Qï (n), sn) Q5 (n) and K(n) using equations 


(4-21) through (4-23) and (4-19). 
3. Set n = 2 and repeat step 2. 


ПА ТЕ кои « m2), скор; п“ =ејљалавеценор та иње cy 


is known. 
5. ІР К(2) « K(1), set n = 3 and repeat step 2. 


6. Continue solving K(n) by increasing n in steps ОЕ 1 
until K(n+l) > K(n). Stop as soon аѕ К(п) < К(п+1); п* = 


current value of n and the optimal solution has been found. 
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The algorithm described in the preceding paragraph was 
used to observe the effect of parameter variations on optimal 
n. In general, the results were similar to those observed in 
the EOQ model. For example, n optimal varied inversely with 
Ar and То; 2 and Ii. 


The sensitivity of n to the various parameters was highly 


but directly with A 


dependent on the selected values of the backorder costs. For 
example, Figure 13 illustrates the relationship between optimal 
n and I, in a situation where the backorder cost is small. 
Figure 14 illustrates this same example with a high time de- 
pendent backorder cost. The difference in the sensitivity is 
obvious. 

Figures 15 and 16 illustrate the relationship of optimal n 
to a change in A): In Figure 15 the backorder costs were low, 


and optimal n was insensitive to a Change in A Та Ерис 


2“ 
example optimal n remained at two, and s* increased as A, 
increased. When the time dependent backorder cost was high 
optimal n varied directly with A, (Figure 16). However, the 
change in optimal n was not nearly as large as in the EOQ model. 
Figures 17 through 20 illustrate the relative magnitude of 


0% ; Q5 , and s* as A, and I, are allowed to vary. The breaks 


in these curves occur at points where optimal n changes. 
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V. RECOMMENDATIONS FOR FURTHER STUDY 


The analysis of the backorders permitted model should be 
completed for all values of the backorder costs T and T in- 
cluding the case where m > 0 and T= 0, Hadley and Whitin 
discuss this case in their single echelon development and 
conclude that if T - 0 then s* is either 0 or infinite. 
Initial investigations indicated that in the two echelon model 
there could exist a finite value of s* » 0 when п = 0. How- 
ever, this investigation was not completed and no conclusions 
were reached. 

It is felt that the analytical argument used throughout 
the paper could be applied to the general backorder case where 
nm and T are both greater than zero. However, the complexity of 
the equations would greatly complicate this development. 

The models discussed in this thesis should be extended to 
more than one activity at each echelon. If it can be assumed 
that all activities within a given echelon order at the same 
time, then it is particularly easy to include multiple activ- 
ities within an echelon. For example, assume that there are 
K activities in echelon one. Then in the EOQ model Ky would 


take the form 


een 
о + EL. 
i-1 11 
If the values of the parameters А, À ir Iii, and C are constant 
Tora Е ey) Kr then. the express он БОК Ky reduces to 
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the trivial case of 


where к, is found in Section II. The expression for K, 
remains unchanged. 

The problem of multiple activities at the top echelon is 
usually not a very interesting one, since each top echelon 
activity is usually responsible for supplying specified activ- 
ities at the next lower echelon. Therefore, unless the model 
provides for lateral resupply actions, the problem of multiple 
activities at the top echelon can be reduced to the sum of a 
series of independent problens, each containing one top 
echelon activity. 

A related problem which should be investigated is the ex- 
tension of the models to more than 2 echelons. In the case 
where each echelon is limited to one activity, this extension 
would not be difficult. Indeed, given the assumption that there 
is an integer ratio between the quantity ordered by successive 
echelons, the model would probably take the form 

m 
KE ии 
ізі. 


where the relationship between the е is 


en en 
йл ен от 
Ф 4-1 
О. = пао ELM = пепо , 
J 2-1 ШІ iz] + T 
and nk 15 an integer for allk, k = L,...,) 1. С Весанзе оние 
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integer property of the п)'5, the solution technique outlined 
in Chapter II could be used to solve recursively for K. 

If the m-echelon problem allows к, activities at the n 
echelon, i = l,...,m, then for large values of к. and m, the 
sheer magnitude of the problem would make its solution ex- 
paenely difficult. Іп fact, it is doubtful that the proce- 
dures suggested in this paper could be utilized for a problem 


of this type. However, if the model is’ small enough, one 


could at least get a "feel" of the model's behavior. 
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VI. SUMMARY 


Three of the deterministic models of Hadley and Whitin 
have been restructured as two echelon models. Equations for 
the average annual variable cost are derived, and standard 
mathematical programming techniques are utilized to find the 
optimal inventory policy. The optimal policy is defined as 
that which minimizes average annual variable cost subject to 
the constraints of the model. 

General solutions are obtained for the EOQ and finite 
production models. Because of the complexity of the cost equa- 
tion the analytical solution of the backorders permitted model 
is limited to the case where m = 0. However, a technique is 
developed which can be used to find the optimal policy in the 
general case, 

A parametric analysis is conducted for each model in order 
to study the behavior of the models under parameter variations. 
The behavior of the curves of n optimal for the EOQ and finite 
production models were identical except that the finite produc- 
tion model enjoyed a significant cost advantage. The behavior 
of the backorder model was very dependent on the magnitude of 
the backorder cost. For example, if the model had a low time 
dependent backorder cost then parameter variations had little 
effect on the optimal behavior of the model. However, if the 
backorder cost was set sufficiently high then the model's 


behavior vas similar to tbat of the EOQ model. Figures 7 


* 
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through 20 illustrate the behavior of the three models under 


varying parameter values. 
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Parameter Values 





units 


I, = 0.75 
C = $100/unit 
70 i 
1 SES $25 
А, = 5200 
n= 5.0 
60 Teac 
A = 100 units/yr. 
50 
40 5 
21 
30 
с Ж 
20 
10 
0.1 0.2 0.3 0.4 0.5 = 


Figure 7 01 р 05 and s* : A, Backorder Model 


with Low Backorder Costs. 
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Parameter Values 


I, = 0.75 
70 С - 51004 Е 
А, = $25 
A, = $200 
п = 5 
п = 100 


100 units/yr 


50 


у 
1) 
ч 
p: 
5 
30 
10 
| ——— d! 
um 092 0.3 0.4 0572 I, 
Figure 18. 07 ; Q5 ала ел” А. Backorder Model with 


High BSchopdor Costs. 
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Parameter Values 


Ii = 0.75 
I, = Ошо 
70 О LOO nlt 
Ај о ита E 
пи = ои Су. 
п = 55.00 
60 т = $15.00 
50 
40 2 
2; 
р 
4:30 
E 
5 
20 * 
m 
* 
10 R 
200 400 600 800 A, 
Figure ac. ОТ ; Q5 Е Е №: I, Backorder Model with Low 


Backorder Costes. 
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